Resonance in Forced Oscillations of an Accretion Disk 
and Kilohertz Quasi— Periodic Oscillations 
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ABSTRACT 

We have performed numerical simulations of a radially perturbed "accretion" torus around 
a black hole or neutron star and find that the torus performs radial and vertical motions at the 
appropriate epicyclic frequencies. We find clear evidence that vertical motions are excited in a 
nonlinear resonance when the applied perturbation is periodic in time. The strongest resonant 
response occurs when the frequency difference of the two oscillations is equal to one-half the 
forcing frequency, precisely as recently observed in the accreting pulsar, SAX J1808. 4-3658, where 
the observed kHz QPO peak separation is half the spin frequency of 401 Hz. 

Subject headings: Stars: neutron — X-rays: binaries — accretion disks — hydrodynamics 



1. Introduction 

Millisecond oscillations in the X-ray light 
curves of systems known to contain neutron stars 
or possibly black holes have been observed for sev- 
eral years with the Rossi X-Ray Timing Explorer 
(see van der Klis (2000) for a review). Recently, it 
has been pointed out that the centroid frequencies 
of the corresponding kilohertz quasi periodic oscil- 
lations (kHz QPOs) are in rational ratios of small 
integers, as 3:2 (Abramowicz & Kluzniak 2001; 
Remillard et al. 2002; McClintock & Remillard 
2003; Abramowicz et al. 2003a, b). This supports 
suggestions that a resonance of some kind is re- 
sponsible for the observed properties (Kluzniak & 
Abramowicz 2001, 2003; Titarchuk 2002). In ad- 
dition, for at least one system in which coherent 
pulsations have been detected, indicating the spin 



frequency (i/^ = 401 Hz in SAX J1808.4-3658), the 
separation in frequency between the kHz peaks 
has been recently reported to be consistent with 
z^s/2 (Wijnands et al. 2003), implying that the 
pulsar is exciting motions in the accretion disk in 
a nonlinear fashion (Kluzniak et al. 2003). 

In this Letter we show that the kHz oscillations 
detected in SAX J1808.4-3658 can be attributed to 
forcing of epicyclic motions in the accretion disk 
by the 2.5 ms pulsar, which induces resonance at 
selected frequencies. The coupling between the 
pulsar and the disk could be due to the magnetic 
field, or to some structure on the surface of the 
star. In other, similar systems, a frequency sep- 
aration equal to the stellar spin frequency is also 
possible. 
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2. Response of a torus to an external per- 
turbation 

In order to study the response of an accretion 
disk to an external perturbation, we make several 
simplifying assumptions. The first is to consider 
not a full accretion disk but a slender torus in hy- 
drostatic equilibrium, orbiting a central body of 
mass M. By virtue of the pressure effects, the ro- 
tation curve is not Keplerian, and for all the cases 
studied here, we consider constant distributions 
of specific angular momentum within the torus, 
which has low mass and is slender in the sense 
that its mass m <C M, and its extension L <^ R, 
where R is the distance separating it from the mass 
M. We thus neglect the self-gravity of the torus, 
which is ellipsoidal in cross section in a meridional 
slice. We additionally assume azimuthal symme- 
try. 

The second assumption is to use a potential for 
the central mass ^kl = M[l — exp{rms/i~)]/fms 
(Kluzniak & Lee 2002), which reproduces two 
features of Einstein's gravity that are relevant 
in our context: (1) there is a marginally sta- 
ble orbit for test masses in circular orbits at 
I'ms = 6M, and (2) the vertical and radial 
epicyclic frequencies (in the equatorial plane) are 
given by C^(r) = (l/47r)^(M/r^) exp(rms/r) and 
K2(r) = (l/47r)2(M/r3)exp(r„,/r)(l - r„,/r), 
respectively, so that their ratio is K{r)/Q{r) = 
(1 — rms/fY^"^ ^ coinciding with the exact solution 
in the Schwarzschild metric in general relativity. 

In hydrostatic equilibrium, pressure and cen- 
trifugal support balance gravity in the radial di- 
rection, and pressure alone balances gravity ver- 
tically. An ideal equation of state for the fluid 
closes the system of equations (we have used an 
adiabatic index 7 = 4/3 throughout). The cen- 
ter of the torus can be conveniently defined as the 
locus of maximum density, where the rotation ve- 
locity is Keplerian, and its position is denoted by 

Clearly a detailed answer to the general prob- 
lem being explored here requires modeling of the 
full accretion disk, whether it is geometrically thin 
or thick. Wc will explore this in future work. For 
the time being, analyzing the behavior of such a 
slender torus can be considered analogous to split- 
ting the disk into thin annuli and investigating 
the properties of one such ring — as in analy- 



sis of radiation induced warping (Pctterson 1977; 
Pringle 1996). The slender torus can be thought of 
as a density enhancement in the accretion disk, the 
properties of which might be transmitted and/or 
amplified into the overall X-ray light curve ob- 
served from the system by an unspecified mecha- 
nism. 

The actual tori are constructed from the speci- 
fied gravitational potential and equation of state in 
standard fashion (see e.g. (Igumcnshchev, Chen & 
Abramowicz 1996; Zurek & Benz 1986)), and their 
dynamical behavior (see below) is followed using 
a two dimensional smooth particle hydrodynam- 
ics (SPH) code (Monaghan 1992; Lee & Ramirez- 
Ruiz 2002) that uses cylindrical coordinates (r, z). 
A cross section (in density) of an equilibrium con- 
figuration thus realized is shown in Fig. 1. We now 
detail the results obtained from using two diflferent 
types of perturbations to such an initial state. 

2.1. Impulsive perturbation 

As a first example, we consider an impul- 
sive perturbation to the torus, analogous to that 
considered recently by Zanotti, RezzoUa & Font 
(2003). (We note that these modes might also be 
relevant for kHz QPOs in black hole systems; sec 
RezzoUa, Yoshida, Maccarone & Zanotti (2003).) 
The equilibrium torus is given a perturbed veloc- 
ity field at f = and allowed to oscillate freely 
thereafter. The velocity field of the perturbation 
is purely radial (in a cylindrical sense). The kick 
is small, so that the induced oscillations in the 
radial direction have an amplitude that is of the 
same order as the extent of the torus itself. The 
torus performs small oscillations in the radial di- 
rection, primarily at a frequency consistent with 
the local epicyclic frequency k,q, as would be ex- 
pected. However, small vertical motions arc also 
induced owing to pressure coupling between both 
modes, and these can be clearly seen in the verti- 
cal oscillations of the center of the torus, occurring 
with greatest power at a frequency that is consis- 
tent with the local vertical epicyclic frequency, Co- 
Both the radial and vertical oscillations arc shown 
in Fig. 2, along with their Fourier transforms. 

It is thus clear that if radial motions of some 
kind are present in the torus, they can induce ver- 
tical oscillations, and both of these occur at the 
epicylic frequencies for test particles in nearly cir- 
cular orbits. 
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2.2. Periodic perturbation 

We now coiisi(ic;i' a situation in which there is 
a periodic perturbation at work, which is exter- 
nal to the disk. The pulsar at the center of the 
gravitational potential well is clearly an example 
of this, and one would expect it to have an effect 
on the disk. The magnetic field of the pulsar, or 
some deformation on its surface, can perturb the 
disk at intervals given by the inverse of the spin 
period, AT = Xjvs (in the case of SAX J1808.4- 
3658, Vs=AS)\ Hz). We have considered then a 
forcing in the radial direction, which manifests it- 
self through a small, radial acceleration, the mag- 
nitude of which is shown in Fig. 3. It is not purely 
sinusoidal but repeats at a fixed interval AT, as 
given above. The adopted profile of this pertur- 
bation deserves some comment: we have used it 
to mimic the passage of a brief (compared with 
the repetition time), but fairly stronger than av- 
erage, disturbance in the accretion disk (e.g., the 
corresponding polar magnetic field of the pulsar, 
sweeping around the disk). Different shapes for 
this pulse will be explored elsewhere. 

The radial forcing obviously induces a radial 
oscillation of the torus. Because of pressure cou- 
pling, a vertical oscillation is again apparent, al- 
beit at a reduced magnitude (compared with the 
radial amplitude). The motion of the center of the 
torus (as defined above) can be Fourier analyzed 
to extract the relevant frequencies. We find that 
the radial and vertical motions occur primarily at 
the local radial and epicylic frequencies, and Co 
respectively, at tq. For a torus orbiting a mass 
M = 1.38Mq, and center at ro = Vl.^'^M = 
24.8 km, Ko = 500 Hz and Co = 700 Hz. 

The question now is: what effect, if any, does 
the repetition time AT of the perturbation have 
on the power of oscillatory motion induced in the 
torus? The radial motion is always driven at a 
relatively high amplitude, simply because the per- 
turbation is itself applied as a radial acceleration. 
We find, however, that the power of the vertical 
motions varies greatly as AT is altered. 

We have kept the initial position of the torus 

fixed at To = 12.25M and performed over a dozen 
simulations, differing only in the value of u^, and 
covering a range 100 < z/s(Hz) < 600. For each 
one of these, we have computed the peak power 
Pm,z in the vertical oscillation at Co- Since one 



could equivalently keep z/s fixed (as is actually the 
case in nature) and vary ro (which would alter the 
difference Co — 'fo)) in Fig- 4 we show Pm,z as a 
fimction of the ratio i^j5/(Co ^ '*o)- The response of 
the torus is clearly greatest when j/^ = 2(Co — ko), 
as was observed in the single instance for SAX 
J1808.4-3658 when two QPO peaks in the kHz 
range were seen. 

We note that there is also a strong response 
when i^s and Co ^ '*o are in a 1:1 or 3:2 correspon- 
dence. The first of these would allow for the possi- 
bility of twin peaks with a separation of 401 Hz in 
SAX J1808.4-3658, while the second would imply 
a separation of 802/3 = 267 Hz. The first option 
would occur at Co=1054 Hz, ko=653 Hz (see also 
Kluzniak et al. (2003)), and ro=9.75M=19.7 km 
for 1.38M0. The second would imply Co=832 Hz 
and «;o=565 Hz at a radius ro=11.13M=22.5 km 
for the same mass. Neither has been observed as 
yet. 

3. Discussion 

The appearance of twin kHZ QPOs in the mil- 
lisecond pulsar SAX J1808. 4-3658, with a separa- 
tion consistent with half the known spin frequency 
of the pulsar, strongly indicates that a nonlinear 
resonance is at work, coupling the spin to vibra- 
tional modes in the disk. Using a simple hydrody- 
namical model, we identify these modes with the 
radial and epicyclic oscillations of fluid elements 
slightly displaced from exact circular orbits. In 
this respect, the model is crucially dependent on 
the effects of strong gravity, to break the degener- 
acy between the orbital and epicyclic frequencies 
present in the Newtonian regime. 

Under the unique assumption that the pulsar 
provides a periodic driving radial force to a slender 
torus in orbit (which we consider as a stand-in 
for a density enhancement in the accretion disk; 
see § 2), we show that the response of the torus 
is greatest when the spin frequency is twice the 
difference between the vertical and radial epicyclic 
frequencies. We also show that there are other 
possibilities for resonant motion, when the above 
numbers are in a 1:1 or 3:2 correspondence. 

For the case of SAX J1808.4-3658, the fre- 
quency ratio 700:500=1.4, and the actual values 
of the frequencies observed would allow us to fix 
the mass of the pulsar at 1.38 solar masses in the 



3 



Schwarzschild metric. However, the actual value 
will be different, as the epicyclic frequencies for a 
neutron star rotating at 401 Hz depart from the 
Schwarzschild values (Kluzniak et al. 2003). 

Two further points deserve comment in the con- 
text of kHZ QPOs in systems with a millisec- 
ond pulsar whose spin frequency is known. First, 
there is an apparent dichotomy in the values of the 
QPO peak separation with respect to the spin fre- 
quency. For the "fast" rotators, like SAX J1808.4- 
3658, the separation is half the spin frequency, 
Az/ = while for "slow" rotators, Uke XTE 

J1807-294 (where i^s = 190 Hz), the separation 
is consistent with the spin frequency, Av = Vg- 
This could be explained within the current frame- 
work, since for fast rotators, the resonant point 
corresponding to Av = is so close to the neu- 
tron star that is it unlikely to be seen (one could 
say that there is no room around the pulsar for 
this mode to occur). For slow rotators, both res- 
onances, at Vg and i^s/2, could in principle be ob- 
served, since they would manifest themselves at 
greater distances from the neutron star. However, 
if the modulations in flux generated in the inner- 
most regions of the accretion disk have a dominant 
effect on the overall X-ray light curve, one would 
then expect to see the strongest signal at frequency 

Vs- 

Second, the twin kHz QPOs then need not nec- 
essarily occur always at the same frequencies, as 
the torus shifts its radial position, but their sep- 
aration should nevertheless remain roughly con- 
stant (and equal to Vg or J^s/2, depending on the 
system) under this excitation mechanism, assum- 
ing mode locking to occur. Further timing obser- 
vations of millisecond pulsars will no doubt shed 
light on these matters. 

More generally, our simulation shows directly 
that a radial perturbation present in the torus 
can excite forced oscillation at (other) eigenfre- 
quencies. This may have implications for mode 
coupling in black hole accretion disks. 
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ments on the original version of the manuscript. 
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Fig. 1. — Logarithmic density contours in a merid- 
ional slice of a slender torus in hydrostatic equi- 
librium, spaced every 0.2 dex, with the lowest one 
at log(p/pniax) = — 1- The center of the torus is 
at ro = 12.25 M. 
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Fig. 2. — Left: Radial (solid line) and vertical 

(dashed line) oscillations performed by the center 
of the torus after an impulsive radial perturbation 
at t = 0. Right: Fourier transforms of the radial 
(solid line) and vertical (dashed line) oscillations 
shown in the left panel. The two curves have been 
rescaled so as to fit on the same graph. 



Fig. 3. — Amplitude of the radial perturbation 
(normalized to the maximum value) applied to the 
torus, with a repetition frequency of z/g = 401 Hz. 
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Fig. 4.-- Peak power in the vertical cpicyclic os- 
cillations induced in the torus by the periodic per- 
turbation of Fig. 3, as a function of Vs/{C,o — i^o)- 
Three peaks at 1.0, 1.5, and 2.0 are clearly visible. 
For this calculation, Vg, the pulsar perturbation 
frequency, varied while Co — '^o was held fixed at 
200 Hz. 



